Quantum Computation with Coherent Spin States and the Close Hadamard Problem 
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We study quantum algorithms inspired by the continuously-parameterized yet finite-dimensional 
Hilbert space of a coherent spin system. We show that the highest-squeezed spin state possible 
can be approximated by a superposition of two states thus transcending the usual model of using 
a single basis state as algorithm input. As a particular example, we show that the close Hadamard 
oracle-decision problem, which is related to the Hadamard codewords of digital communications 
theory, can be solved quantitatively more efficiently using this computational model than by any 
known classical algorithm. 
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I. INTRODUCTION 

Ever since the remarkable discovery that quantum me- 
chanical systems can in principle be used for computa- 
tional purposes by pioneers such as Benioff [ij , Feynman 
, and Deutsch Q , researchers have considered the feasi- 
bilities of implementing such computations in a rich vari- 
ety of physical models. These models include phenomena 
such as light, electrons, atomic nuclei each of which have 
degrees of freedom such as position, path or spin into 
which quantum information can be encoded, processed 
and measured. The measurement spectrum can be dis- 
crete, for example measuring spin up or down or position 
being left or right, or continuous, for example measuring 
the angle of the spin axis or where on a line the particle 
is located. 

This distinction between whether the spectrum of mea- 
surement, or preparation, is discrete or continuous is at 
the heart of the difference between discrete variable and 
continuous variable versions of quantum information. A 
system is considered to be continuous variable when ei- 
ther the spectrum of preparation or the spectrum of mea- 
surement is continuous. 

The processing of quantum information in discrete 
variable systems has the advantage of being amenable 
to the theory of quantum error correction. Discrete 
variable systems, where experimentalists have success- 
fully demonstrated the potential computational power of 
quantum systems, include trapped ions liquid nuclear 
magnetic resonance Q and quantum dots @ . 

Processing with continuous variable systems has the 
advantage of being amenable to many physical prepara- 
tion and measurement procedures that feature continu- 
ous tunability. Continuous variable quantum informa- 
tion is the less studied of the two types yet offers new 
insights that can engender creative thinking about new 
quantum algorithms. Given the significant progress that 
has been made in recent years in understanding and con- 
trolling quantum optics systems [7n-15j. there seems to 
be a potential for further progress in realizing quantum 
computation using continuous variables. 



Continuous variable studies are typically linked to har- 
monic oscillators because quantum optics has powerful 
tools to work prepare, process and measure optical field 
modes ^1^, which are analogous to harmonic oscilla- 
tors. Quantum optics also offers an easy-to- implement 
Fourier transform, which maps canonical position states 
to canonical momentum states. The duality between po- 
sition and momentum states is evident in that the magni- 
tude of the overlap between any position and any momen- 
tum state is constant regardless of the choice of position 
and momentum states. 

This duality is only exact for continuous position and 
momentum over an unbounded spectrum, which is an 
unattainable ideal. A well-know example of this dual- 
ity is the time-bandwidth uncertainty relation of Fourier 
transforms. This inherent uncertainty relation limits the 
performance of quantum algorithms implemented over 
qudits encoded into position states of a single harmonic 
oscillator [iSi] . Encoding over sufficiently many har- 
monic oscillators [Toj ameliorates this limitation some- 
what. 

Continuous variable quantum information processing 
may be thought of as a model of quantum computation. 
The circuit model was the first such model [2(|, and the 
circuit and one-way quantum models are examples of dif- 
ferent models [2l|, which are distinguished by different 
processing methods. The continuous variable model of 
quantum computation has inspired a continuous variable 
Deutsch- Jozsa algorithm 0, [13, Es]. Here we explore an- 
other continuous variable model of quantum computation 
based on coherent spin states and show that it can inspire 
a new algorithm. 

Although fixed total spin states span a finite- 
dimensional Hilbert space, continuous encoding is pos- 
sible using the continuously-parameterized, squeezed 
spin coherent states [22], which are analogous to the 
squeezed coherent states of the harmonic oscillator (23l - 
|25| . The finite-dimensional Hilbert space has the ef- 
fect tha t sq ueezing of coherent spin states is Heisenberg- 
limited [22| , unlike the coherent states of the harmonic 
oscillator where the allowed amount of squeezing is un- 
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bounded. We consider encoding quantum information 
into the highest-squeezed spin state that can be achieved. 

We demonstrate that this optimally squeezed state 
may be approximated well by a superposition of two dis- 
crete states, thus idealizing the computation model be- 
yond encoding into squeezed spin states while keeping 
the spin gates. This approach allows us to discover a 
new algorithm that can be processed using the circuit 
model of quantum computation. Our investigation of a 
continuous variable spin model of quantum computation 
has thus inspired us to find a new quantum algorithm. 

Our paper is presented as follows. In Sec. II, we 
introduce an oracle decision problem parameterized by 
N = 2"-bit strings. We refer to this new oracle de- 
cision problem as the close Hadamard problem, which 
is related to the digital coding techniques employed in 
classical communications [1^ . The close Hadamard 
problem is inspired by the continuous variable opera- 
tors and displacement and squeezing tools used in the 
continuously-parameterized, infinite-dimensional Hilbert 
space case adapted to the continuously-parameterized, 
finite-dimensional Hilbert space case. 

In Sec. HI, we introduce the spin-system model. We 
discuss the preparation of coherent spin states and how 
to use Q-functions for spin systems, parameterized by 
the spherical angles, as a visualization aid. We intro- 
duce the concept of spin squeezing and show how squeez- 
ing changes the amplitude distribution of the individual 
spin states. We show, that for a particular coherent spin 
state, the limiting squeezed state is asymptotically ap- 
proximated by a symmetric superposition of two discrete 
states with constant error independent of the size of the 
Hilbert space. We use this superposition as the algorithm 
input state. 

In Sec. IV, we demonstrate that our quantum algo- 
rithm efficiently solves the restricted close Hadamard 
problem with certainty in a single oracle query and also 
solves the unrestricted close Hadamard problem with ar- 
bitrarily small error in a constant number of queries. 
We also show that any known classical algorithm re- 
quires f2(n) queries. In the restricted case in particular, 
this speed-up is the result of the cancellation of bit errors 
of certain patterns and results from using a superposition 
of two states as algorithm input. 

In Sec. V, we discuss generalization of the computa- 
tional model by showing that if the Hadamard operation 
is replaced by the discrete Fourier transformation, the 
oracle decision problem changes and is related to what 
can be thought of as fractional bits. We conclude in Sec. 
VI that this model of quantum computation can be used 
to inspire the efficient solution of additional problems. 



II. ORACLE DECISION PROBLEMS AND THE 
CLOSE HADAMARD PROBLEM 

Quantum algorithms for the efficient solution of or- 
acle decision problems are of historical importance in 
quantum information 0, . The Deutsch- Jozsa prob- 
lem in particular has been studied in both discrete and 
continuous variable settings 0, [13 [3 HI]- We are 
inspired by the continuous variable quantum algorithm 
used to solve the Deutsch- Jozsa problem, where logi- 
cal states are a continuously-parameterized in a infinite- 
dimensional Hilbert space. Here we explore this continu- 
ous variable quantum algorithm where the logical states 
are continuously-parameterized in the /inite-dimensional 
Hilbert space of a spin system. 

This continuously-parameterized spin model of quan- 
tum computation naturally yields to a symmetric super- 
position of two basis states as the logical state. Returning 
to a discrete representation inspires us to discover a new 
oracle decision problem, which we refer to as the close 
Hadamard problem. 



A. Oracle Decision Problems 

Oracle decision problems are related to oracle identi- 
fication problems, which are usually presented in terms 
of the problem of identifying a unique function /. The 
function / maps N = 2"-bit strings to a single bit 

/: {0,1}" .^{0,1}. (1) 

Any Boolean function on n bits can also be represented 
by a string of = 2" bits, in which the i^^ bit Zi is the 
value of the function on the i"^ bit string, taken in lexico- 
graphical order. The challenge of the oracle identification 
problem is to identify a unique iV-bit string from a set of 
size 2^ by making the fewest queries to an oracle. 

Oracle decision problems are simpler than oracle iden- 
tification problems because the function or string does 
not have to be identified explicitly. Rather the problem 
is to identify which of two mutually disjoint sets contains 
the string. For our analysis, the oracle decision problem 
is defined as follows. 

Definition 1. An oracle decision problem is specified by 
two non-empty, disjoint subsets A, B C {0,1}^. Given 
a string z G AU B ^ C , the oracle- decision problem is 
to determine whether z d A or z d B with the fewest 
queries to the oracle possible. 



B. The Close Hadamard Problem 

In this subsection, we specify the particular sets A, B 
and C required by Definition 1 for the close Hadamard 
problem. We refer to this decision problem as close be- 
cause we are interested in strings that are close in the 
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sense of Hamming distance to the A'' 
referred to as Hadamard codewords [H, |29| . 

The problem of discriminating between codewords re- 
ceived after transmission over a noisy channel is well- 
known in classical digital coding theory employing linear 
block codes [sO] • Linear block codes are characterized by 
the triplet [N,k,t], where N is the total length of the 
codeword; k < N is the amount of information coded, 
and t — 1 is the number of errors that the code can cor- 
rect. 

The Hadamard code is a linear block code with N — 
2", A; = n + 1, and t = N/i ~ 1. The Hadamard code 
has a poor information rate fc/A^, but it has excellent 
error-correcting capability. Because of this latter feature, 
the [32, 6, 7] Hadamard code was used to encode picture 
information on Mariner space craft missions [29i] . 

For N — 2", the matrix comprising Hadamard code- 
words is 



2"-bit strings codeword wf^. We define the set of codewords with 



W(^) = log, 



(-1) 



7VH 



(2) 



where H**" is the familiar Hadamard matrix of quantum 
computation [3l|. The expression log(_i)[a;] is the log- 
arithm of X to the base —1, and, if a; = then 
y = log(-i)N. 

For counting purposes we define the set Zjv = 
{0, 1, . . . , iV - 1}. For 2 G Zat, the j*'' Hadamard code- 
word corresponds to the f"^ row of the matrix W^^' and 
is expressed as w\^^ . For example 





10 1 

11 

110 



(3) 



and W.. 



(4) 



0110. Note that all Hadamard codewords are 
balanced with the exception of wl^\ which is constant. 
Also note that all N Hadamard codewords are separated 
from each other by Hamming distance 



errors specified with respect to any particular codeword 
through the use of an error syndrome, which represents 
all the possible ways an error of d bits can occur. 
The error syndrome for d unrestricted errors is 



f/, = {ze{0,i}^||2| = d} 



(5) 



The set of codewords having d unrestricted errors with 
respect to the j*^ codeword is 



(6) 



and the set of all correctable codewords with zero to A^/16 
unrestricted errors is 



(7) 



We proceed in a similar manner with the definition of the 
Hadamard codewords having restricted errors, which are 



a subset of S 



2. Codewords with restricted errors 

The codewords having restricted errors are the strings 
with Hamming distance d from Hadamard codeword 
Wj'^\ but the errors are restricted to the N/2 specific 

bit positions where the codeword wj^J^ contains a one. 
The error syndrome for d restricted errors is 



Rd^{ze {0, 1}^ \\z\^d and z ^ wj^^\} 



(8) 



We say that a vector a E {0, 1}^ is dominated by a vector 
b e {0, 1}^, denoted a ^ 6, if whenever = 1 then also 
bi ~ 1. The set of codewords having d restricted errors 
with respect to the j"^ codeword is 



• w. 



z e 



Rd} 



(9) 



N/2. 



(4) 



We present two examples of the sets given by Eq. ([9]). 
For the N = 8 case where there is a single restricted 

An arbitrary string z G {0, 1}^ having Hamming dis- error, we have Wy*^ = 01101001, and, for the particular 



tance d{ z,W, 



' ^z, Wj^''^ < N/A from any Hadamard codeword 
is said to be within the error-correcting capability of 
the Hadamard code [27j . 

In our analysis, we introduce Hadamard codewords 
with two types of bit errors: unrestricted errors and re- 
stricted errors. Unrestricted bit errors can occur at any 
of the N bit positions, whereas restricted errors are lim- 
ited to N/2 specific bit positions. 



1. Codewords with unrestricted errors 

The codewords having unrestricted errors are the 
strings having Hamming distance d from any Hadamard 



codeword w]®^ = 00001111, 
Efl ={01001111,00101111,00000111,00001110}. (10) 



Inspection of the set given in Eq. (flO)) reveals the error 

(s) 

alignment with the bit positions where Wj =1. 

For iV = 8 where there are two restricted errors, the er- 
rors may occur at any two of four possible possible bit po- 
sitions represented as {a, b, c, d}, for which there are the 
(2) = 6 distinct bit error pairings {ab,ac,ad,bc,bd,cd}. 
The codewords with two errors in this case are 

Efl ={00000110, 00100111, 00101110, 01000111, 



01001110,01101111}. 



(11) 
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For the general case, the set having m-tuple restricted 

errors has size S^-^^ = (^^) . 

The set of all correctable codewords with zero to N / 4 
restricted errors with respect to the j*'' Hadamard code- 
word is 



The size of this set is exponential in N since 



(12) 



N/i-1 
m=0 



2^ 



(13) 



We now define two variations of the close Hadamard 
problem in terms of Definition [TJ 

Problem 1. Given the set of codewords A — sj^2_i) 
which contains strings that are close (in the restricted 
sense) to the Hadamard codeword wjyl_-^ and the set 

of codewords B — S^,^'' , which contains strings that are 
close (in the restricted sense) to any other Hadamard 
codeword wjf^^ with k G ^nd a string z randomly 

selected with uniform distribution ^ such that z C = 
AU B, the restricted close Hadamard problem is 
to determine if z £ A or z & B with the fewest oracle 
queries. 

In our formulation of Problem [U we have made a tech- 
nical assumption by setting j = N/2 — 1 in our definition 
of set A. We could have selected any other j G as 
long as we excluded the selection from the definition of 
set B. We make this assumption because our quantum 
algorithm requires the measurement of some qubit. We 
have arbitrarily, and without loss of generality, set it to 
the qubit that corresponds to j — N/2 — 1. The same 
assumption is made in our formulation of Problem [2j 

Problem 2. Given the set of codewords A = S^2_ij 
which contains strings that are close (in the unrestricted 

and the set 



sense) to the Hadamard codeword 



(AT) 

N/2-1 



of codewords B = , which contains strings that are 
close (in the unrestricted sense) to any other Hadamard 
codeword wjf^"^ with k G '^n/2-i ^^d. a string z randomly 
selected with uniform distribution /i sucii that z G^ C = 
AU B. The unrestricted close Hadamard problem 
is to determine if z £ A or z £ B with the fewest oracle 
queries. 



C. Quantum Algorithm for Oracle Decision 
Problems 



The quantum circuit represented in Fig. [T] solves the 
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FIG. 1: Single-mode quantum circuit for oracle decision prob- 
lems fil. 



Deutsch-Jozsa oracle decision problem employing logical 
states encoded in the infinite-dimensional Hilbert space 
of the harmonic oscillator [17, 18]. In this case, 4'(x) is a 
square- integrable function of the continuous position x. 

A key aspect of the approach is the physical acces- 
sibility of the harmonic oscillator ground state and the 
availability of linear and quadratic operators that enable 
us to prepare the logical input state |5'o). Additionally, 
the operators R and identified in Fig. [T] are the eas- 
ily implementable continuous Fourier transform and its 
inverse. The oracle information is encoded in the logical 
state by the unitary operator [l3l 



/ i-iy^ 

(-1)^= 



V 







\ 





(14) 



where the Zi are the bits of the unknown string given in 
oracle decision problem Definition 1. 

Mapping this algorithm to the finite-dimensional, 
continuously-parameterized coherent spin system must 
first deal with the step that takes the ground state of 
the spin system to the logical input state jvPo)- This first 
step also employs physically accessible, linear spin rota- 
tion and quadratic spin squeezing. 



III. THE SPIN SYSTEM MODEL 

The method of generalized coherent states has been 
successfully used to describe a number of diverse physi- 
cal phenomena including quantum optics, atom-light in- 
teractions, and superfluidity [23j . Here we make use of 
coherent spin states [13, Hal in creating an alternative 
model of continuous variable quantum computation. Just 
as squeezing is beneficial in continuous variable quantum 
computing using the harmonic oscillator, we make use of 
spin squeezing here [1^. We use the optimally squeezed 
spin state [1^ as input to our algorithm and show that it 
can be approximated by a superposition of two discrete 
states with constant error independent of the size of the 
Hilbert Space. 



5 



A. Coherent Spin States 

Our spin system is a collection of 2S elementary 1/2 
spins. 25 is thus an odd integer, and we choose 25* + 
1 = iV so that N = 2"-bit strings may be naturally 
represented. We refer this as an S'-spin system [2^ . 

The system dynamics are determined by the Hamilto- 
nian, which is expressed as a polynomial of su(2) alge- 
braic elements. These algebraic elements are Pauli spin 
operators in the spin-1/2 single-particle case. For higher 
even dimensions, we use notation similar to j2^ with op- 
erators Si, Sj and Sfc and i, j, fc denoting the components 
of any three orthogonal directions, such that 



S S 



and 



AS?AS|>^(S.) , 



(15) 



(16) 



and cyclic permutations. 

The spin system is oriented in the usual way. With 

me {-s,-s + l,-s + 2,...,s}, (17) 

the spin kets \m) ^ are eigenstates of Sz and satisfying 

Sz = m |m)^ , (18) 

and 

S2|m)^ = s(s + l)|m)^, (19) 

where = + + S^. The ladder operators are S± = 

Sx ± iSj, , and the action of the lowering operator on the 
ground state is 



S-l-s) =0. 



(20) 



We use the discrete spin states to construct the continu- 
ously parameterized coherent spin states. 

The harmonic-oscillator coherent states are transla- 
tions of the oscillator ground state [l6|. Analogously, 
the coherent spin states are rotations of the spin-system 
ground state 22-24]. Individual spin states are often 
referred to in the literature [l^, H^] as Dicke states 
analogous to photon number states, and the coherent 
spin states are referred to as Bloch states analogous to 
Glauber states. 

The coherent spin state, 16*, 0) with 6*, e M, is [22j |: 



1^, 



Re 



= 1 + tan^ - 



2s 



E 

fc=0 



e"* tan ■ 



(21) 




0.1 




(b) 



64 



FIG. 2: 

for s = 



(a) Spherical Q-function of the state given by Eq. (I22p 
and (b) Plot of the respective Dicke-state proba- 



bility distribution. 



The coherent spin state of most interest is 



|7r/2,0), = 2-£ 



k=0 



(22) 



which has a Dicke-state amplitude spectrum whose 
squared magnitude is the binomial probability distribu- 
tion with p ^ q = 1/2 shown in Fig, ^jp. 

Quasi-probability distributions [23| are a useful means 
of visualizing spin states. We choose to use Q- 
functions [23| as coherent-state representations. The 
spherical plots of these distributions provide good intu- 
ition as to the orientation and the isotropic or anisotropic 
distribution of uncertainties, but they are not used to ac- 
tually calculate uncertainties, which is performed using 
Eq. ([HI). 

For the arbitrary coherent spin state represented as 
\^\ ~ 'Ylik=Q'^k we express the spherical Q-function 
llal as 



2s 



2s 
fc=0 



sin(6l/2)''cos(6'/2) 



2s- 



akc ^. 



(23) 



In Fig. [21 we plot the spherical Q-function and the prob- 
ability distribution of the discrete spin state for the state 
given in Eq. fF2\ . Note that this coherent spin state ap- 
pears as an 'equatorial' state with isotropic uncertainty 
distribution when represented this way. 



B. Squeezed Spin States 

Coherent spin states can be squeezed [22] . Whereas the 
Glauber states can be squeezed to an arbitrary degree, 
spin states can only be squeezed to the Heisenberg limit 
of 1/2 . We wish to exploit the squeezed state with 
the minimal achievable variance in our algorithm. In 
the following, we formulate expressions for this optimally 
squeezed spin state and show that it can be approximated 
well by a superposition of two spin states. 
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FIG. 3: (a) Spherical Q-function of the squeezed state given 
by Eq. ((ST]) for s = ^, and (b) Plot of the respective Dicke- 
state probability distribution. 



We employ two-axis counter-twisting [22| to define the 
squeezing operator 

=e'^^.giA.(s^-SJ)^ (24) 

where fj, is the squeezing parameter [22j . The rotation op- 
erator e'?'^" orients the resulting anisotropic uncertainty 
distribution in the y, z directions. 
Applying the operator to 



l*) = k/2,o). 



(25) 



allows us to reduce the variance AS^ at the expense of 
enchancing the variance AS^ . The reduced variance may 
be expressed as 



V- 



(26) 



since the first moment (Sz) = 0. In Fig. [3^, we plot the 
quasi-probability distribution of a squeezed spin state. 
The reduced variance of the squeezed state in the z direc- 
tion and increased variance in the y direction is evident. 

The minimum value of the reduced variance V- asymp- 
totically approaches 1/2 with increasing s [l^. We re- 
fer to the optimal value of the squeezing parameter at 
this minimum as /iopt- For /i > /iopt, the distribution 
variance increases and the distribution quasi-probability 
distribution becomes becomes skewed [22| ■ It can shown 
be that fiopt — >■ 7 as s — > 00. This limit is understand- 
able since the variance of a binomial distribution with 
p = q — 1/2 is N/A, and squeezing simply has the effect 
of removing the distribution variance of the dependency 
on iV = 2s + 1. 

We express the optimally-squeezed spin state as 



<i>W\ = |V2,0,A*opt), 



(27) 



with I*) defined in Eq. (gHI)- In Fig. [3Jd, we plot 
the Dicke-state probability distribution of the optimally 
squeezed state. It is evident that this state approximates 



the superposition of two spin states. We wish to provide 
a bound on how well approximated the squeezed state is 
by a two-component superposition. 

Analysis of the variance of the squeezed state's proba- 
bility distribution is facilitated using the qudit represen- 
tation rather than the spin sate representation. We thus 
represent this iV-dimensional squeezed state in terms of 
the qudits |i) as 



N-l 
i=0 



(28) 



The probability distribution associated with |$(^)) may 
be represented as the set 



\aN- 



(29) 



with individual probabilities V^^^ = We note that 

the squeezed state is symmetric about the centre, and 
thus the two central states have 



,(JV) 
{N/2-1) 



(N) 
{N/2: 



- piN) 
^ c ■> 



(30) 



and thereby form the principle components of the prob- 
ability distribution of optimally squeezed states. 

For s > |, the expression for the reduced variance 
given by Eq. ()26|) requires solving eigenvalue problems 
of degree greater than eight and is no longer analytic, 
and we must resort to numerical analysis. For s = 3/2 
and = 4 the expression for the reduced variance given 
by Eq. (f26|) is analytic, and /iopt — We can thus 

represent the optimal squeezed state as 



c&W \ = e'M 0|0) + -L|l) + -L|2) + 0|3) ) , (31) 



where e^^ is a global phase picked up by the action of . 
The associated probability distribution is 



75(4) ^ {0,1/2,1/2,0}. 



(32) 



For this case we achieve what we refer to as 'perfect' 
squeezing, where 'perfect' means that the two central 
components have probability equal to a half, and the 
probability of the other two components is zero. 

However, this four component distribution has a vari- 
ance of only a quarter, where the distribution variance is 
expressed as 



N-l 



Var 



■2 I |2 

I \ai\ - 



(33) 



Indeed for aU A^ if V?^^ = 1/2 then. 



Var 



^{N) 



J-{{Nl2-lf + {N/2f)^\{N~lf 



(34) 
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FIG. 5: Histograms of the probability distribution p'^' for 
(a) s = 7/2, (b) s = 15/2, (c) s = 31/2 and (d) s = 63/2 with 
a logarithmic scale for the ordinate. The bounding distribu- 
tion given by Eqs. (|36|l and (|37|) is overlayed on each of the 
histograms. 



FIG. 4: (a) Calculated value of the reduced variance of the 
probability distribution given by Eq. (1261) approaches 1/2 with 
increasing s as predicted 22]. (b) Calculated value of the 
probability of the two central components given by Eq. (|30p 
is bounded by the constant given by Eq. (|37|l . 



Since the distribution variance approaches 1/2 as N = 
2s + 1 approaches infinity, perfect squeezing in the sense 
we have defined is not possible. We use the variance 
equals 1/2 as a means to bound V^'^^ defined in Eq. (l30l) . 
In Fig. |4^, we plot the calculated values of the reduced 
variance given by Eqs. (|26|) and (p3|) as a function of s 
from s = 3/2 to s = 1023/2, where we observe that the 
variance approaches V — 1/2 as predicted. In order to 
bound the Hmiting value of the two central components 
Vc^\ we bound the 'tails' of the probability distribution 

In Fig. O we plot histograms calculated from the 
squeezed distribution, 'P^'^\ for several values of s, where 
we have scaled the ordinate to reveal the structure of the 
tail components. We see that the components immedi- 
ately adjacent to the central components have 



' N/2+1 ~ ' N/2~2 



0, 



(35) 



and further outlying terms tail off in an exponential-like 
fashion. We thus introduce the following bounding prob- 
ability distribution 



[N] 



e 2e 1 

■■'«'3'y'«'2 

2e e 

.,o,-,-,o. 



(36) 



..,0 



in order to calculate a bound on the two central compo- 
nents of the distribution. 



Solving Var 



= 1/2 for e gives the probability of 



the two central components 



0.484. 



(37) 



In Fig. we plot the calculated values of Vc^\ where 
we note that it goes from 1/2 at s = 3/2 and asymptot- 
ically approaches the constant bounded from below by 
Eq. p7p . The bounding distribution is also overlayed on 
the histograms presented in Fig. [S] 

Since greater than 98% of the probability is manifest 
in the two central components, we approximate the op- 
timally squeezed input state by the superposition of two 
states 



V2 ( 



(38) 



and we can leave the spin system model behind because 
the ideal state given in Eq. (p8| can be used as input 
on any quantum computer. We emphasize that the idea 
to use it as input to the algorithm give in Fig. [T] has 
been inspired through the analysis of the continuously 
parameterized finite-dimensional Hilbert space of the s- 
spin system. 



IV. EFFICIENT ALGORITHM FOR THE 
CLOSE HADAMARD PROBLEM 

An important principle in quantum information pro- 
cessing is the solving of problems with increased efficiency 
compared to classical information processing. Efficiency 
can be measured in terms of a problem's query complex- 
ity, and in this section we present two claims quanti- 
fying the query complexity required to solve the close 
Hadamard problem in the quantum setting. 
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For a quantum algorithm employing the quantum cir- 
cuit in Fig. [T]with R = Rt = H®", we state the number 
of oracle queries required to solve the restricted and the 
unrestricted close Hadamard problem defined in Prob- 
lem [T] and Problem [2] respectively. We prove these claims 
through the introduction and proof of three Lemmas. We 
follow this with a discussion of the query complexity of 
classical algorithms. 

With the idealized input state simplified to 



(39) 



where we have suppressed the normalization factor in 



Eq. ((38)) . the action of the algorithm on the input state 
is expressed as 



l^-,) =H®"U,H 



(40) 



The A''-bit string z represents the function /. We show 
that the algorithm efficiently solves both versions of the 
close Hadamard problem. 

Claim 1. The restricted close Hadamard problem can be 
solved with certainty in a single oracle query using the 
quantum circuit given in Fig. [7] and the idealized input 
state given by Eq. i39\). 

Claim 2. The un-restricted close Hadamard problem can 
be solved with arbitrarily small error probability in a con- 
stant number of oracle queries using the quantum cir- 
cuit given in Fig. [I] and the idealized input state given by 
Eq. M 



Algorithm Response to the Hadamard 
Codewords 



A key feature of the input state is that it is a symmet- 
ric superposition of two basis states. When Hadamard 
codewords are encoded into the oracle, the action of the 
algorithm preserves the symmetric superposition. This 
preservation is demonstrated in Fig. [51 We thus prove 
the following Lemma is true. 

Lemma 1. Given the input |5'o) = |^)^ -I- | — to 
the circuit shown in Fig. Q] and the oracle encoded with 
one of the Hadamard codewords z = W'j^'* for N = 
2s -I- 1 and < j < output state is another super- 

position of spin states l^a) = |i + j)^ + ~ j)^- 



Proof. In order to simplify notation in the following, we 
suppress the superscript N in W^^^ . With W defined as 
the set of Hadamard codewords given by Eq. ([2|), the pair 
(W, ©) forms a group under addition modulo two [27'|. In 
particular, the identity element is Wq, and each element 
is its own inverse since 



W, ® Wj = Wo. 



(41) 



H'3> = |l/2>7/2 + |-l/2>7/2 (a) 



1 2 3 4 5 6 7 



|i> 



H'3>=|3/2>7/,+|-3/2>,/2 (b) 



H'3> = |5/2>7/2 + |-5/2>v,2 (c) 



1 2 3 4 5 6 7 



|i> 



1 2 3 4 5 6 7 



H'3>=|7/2)7,2+|-7/2>,/2 (d) 



1 2 3 4 5 6 7 



|i> 



FIG. 6: For the example oi s — 7/2, A'" = 8, the output state 
given by Eq. (|43p remains a symmetric superposition of two 
states for the Hadamard codewords: (a) Wg*', (b)Wj*', (c) 
, and (d) W^^\ 



It follows from the group property that the addition 
of any two codewords is another codeword. For the 
Hadamard codeword pairs Wj and W^-i-j, it can be 
readily shown that 



Wj ® Wn-1-j = Wn-1 



(42) 



for all j £ In- 

With some algebraic manipulation, the state j^Pa) may 
be expressed as 



y=0 \x=0 / 



(43) 



We use the qudit representation \y) rather than the spin 
state representation \m)s. We translate back to spin state 
representation as the last step. 
The symbol 

=(-l)(^®^«+"/2). + (_i)(^®i^iv/2-(„+i))^^ (44) 

where z — Wj is the string encoded in the oracle, and 
the symbol x represents the a;**' bit of the TV-bit strings. 
Note that the sums y + N/2 and 7V/2- (y-l-1) in Eq. ((ii)) 
are modulo N sums. 

The Hadamard codewords are balanced with the ex- 
ception of Wo, which is constant. For j ^ k, this allows 
us to write 



iV-l 



and for j 



x=0 



N-1 



(45) 



(46) 



x=0 
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where we have used the group inverse relation given in 
Eq. (mj. Using this result, we see that a non-zero sum of 
ax,y in Eq. (j43|) occurs exactly twice when y = {N/2+j) 
mod N and when y — {N/2 — 1 ~ j) mod N. 

In the qudit representation, the output state is thus 
expressed 



1*3) = |iV/2-l-j) + |iV/2 + j) 



(47) 



where qudit indices are understood to be modulo N. For 
< j < N/2, this translates back to the spin basis as 



1*3) ^ 1-1/2 -j->^ + |l/2 + j)^ 
thus completing the proof of Lemma 1. 



□ 



We now show that the superposition of two states is 
also preserved for codewords with restricted errors. 

Lemma 2. Given the input |5'o) = |^)^ + 1^-^)^, to 
the circuit shown in Fig. Q] and the oracle encoded with 
z G sj^"* given by Eq. which is a codeword having 

restricted errors, and with < j < ^, the output state 
is another superposition of spin states [^'3) = | ^ + i) + 



Proof. Observe that the form of Eq. p4)) allows for can- 
celling of errors. Under certain conditions if an error 
occurs at bit position a:, the effect on the left-hand side 
of the plus sign is cancelled by the opposite effect on the 
right-hand side. Consider the x*^ bit error in Eq. (|44)). 
where we have cancellation 

^(_i)(^®^«+«/2)^ _l_ (^_2)(^®^"/^"<"+^')x. (49) 
This identity implies that 

1 = {Wy+N/2 ® WN/2-{y+l))^ 

= {Wn-i)., (50) 



where we have used the result expressed in Eq. (j42l) . This 
is exactly the same as the requirement to be a member 
of set S^^-* defined in Eq. (|12p . Under this condition, the 
result of Lemma 1 holds and I ^E* 3) ~ I5 + + ^ j),- 



□ 



We now show that the perfect superposition of two 
states is no longer preserved for codewords having unre- 
stricted errors. The effect of errors that are not of the 
restricted type is to degrade the superposition by dis- 
tributing amplitude evenly across all other states. How- 
ever as long as the number of these errors is less than 
iV/16, it is still possible to efficiently identify the desired 
state. 

Lemma 3. Given the input \'^q) = \^) ^ + \~^) ^ to 
the circuit shown in Fig. Q] and the oracle encoded with 
z g '^j^'* given by Ec^. dJ), which is a codeword having 
less than N/16 unrestricted errors, the desired state can 
be identified with success probability of at least . 



Proof. Let Vj G ^fi^ be a Hadamard codeword with a 
single unrestricted error. We have already shown that 
if the error is of the restricted type, two-component su- 
perpositions are preserved. If the the error is not a re- 
stricted error, we have no error cancellation, so the single 
bit-error breaks the balanced and constant sums defined 
by Eqs. pS)) and respectively. For j ^ k this gives 



N-l 

2^ ^(_l)(w,eyfe). 



a:=0 



(48) and for j — k, 



N-l 

7V-2 = ^(-l)(^^-®^^)=^. 



(51) 



(52) 



As the input state is a two-component superposition, the 
above sums result in all the amplitudes of the output 
state either acquiring or losing an amount of amplitude 
proportional to four — two from the amount in the bal- 
anced or constant sums given in Eqs. ((5T|) and (|52p and 
two from the effect of there being two components in the 
input state. 

Thus, we express the output state for the worst case 
of a single unrestricted error as 



1 



V2N 



k=s~l/2 

E ± 

k = -s-l/2 



(53) 



For the worst case of I unrestricted errors, where no errors 
are of the restricted type, the principle components have 
amplitude 



If Al\ 



and the amplitude of the next largest component is 



(54) 
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V2N' 



(55) 



The amplitude reduction of the principle components by 
an amount directly proportional to the number of errors 
results from the constant sum given by Eq. (j52p being 
reduced by double the number of errors. However, the 
balanced sums are variable since errors can cancel. The 
worst case occurs when the errors are 'in phase' result- 
ing in the amplitude of the next-largest component being 
proportional to the number of errors. The effect of code- 
words with unrestricted errors on the input superposition 
is presented in Fig. [T] 

The probability of the two central components is 



1 - 



8/ 

N 



16P 



(56) 
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(d) 
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FIG. 7: For A'^ = 64, the effect of I unrestricted errors, where 
none of the / errors are of the restricted type, on tlie proba- 
bility of the central components is demonstrated for (a) I = 0, 
(b) 1 = 2, (c) I = 4, and (d) I = 6. Fig. 7(c) corresponds to 
the case that I = 7V/16. 



The equality holds for the worst case where none of the 
errors are of the restricted type. If I = N/16, then \a\^ > 
i (jg). The amplitudes of the two central states can be 

combined into a single state with amplitude V2a by an 
appropriate unitary operation. Since I is less than A^/16, 
the desired state can be identified with probability 2\a\ , 
which is at least f^. □ 
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B. Solution of the Close Hadamard Problem 

We now use results of Lemmas 1, 2 and 3 to show 
that Claims [1] and [5] are true. 



If j = N/2 - 1, then Pr[s] = 1 and z e A, and if 
j ^ N/2 - 1, then Pr[s] = and z g B. Thus, 
the restricted close Hadamard problem is solved with 
certainty in a single oracle query. □ 



Proof of Claim I 



By Lemmas 1 and 2, I'^s) = 



z e ^ 



A(JV) 



-4 



for 



By Lemma 3, the effect of including I unre- 
stricted bit errors on the output state may be expressed 
as 



+ J) + 



7E± 

A 



(60) 



where the symbols a and /3 are given by Eqs. (|54p 
and (|55|) respectively, and I7I < |/3|. Note that k + A = 
N — 2, so that all N possible states are accounted for, 
but the specific value of 7, k and A are dependent on N 
and I. The outcome probabilities of measuring the state 
Is), are thus 



Prfsl 



1 



+ i 



1 



+ i 



(61) 



If j = N/2 - 1, then Pr[s] > la] , and if j ^ N/2 - 1, 
then Pr[s] < |/3|^. Assuming that the number of unre- 
stricted errors I is less than A^/16, then an error of 0(e~') 
can be achieved by making 0{q) repetitions of the algo- 
rithm 17|. Thus, the unrestricted close Hadamard prob- 
lem is solved with arbitrarily small error probability in a 
constant number of queries. □ 



Proof of Claim [T} 

By Lemmas 1 and 2, l^fs) = \ \+ + \-\ - j) ^ for 
z e Immediately prior to the measurement step. 



we require a unitary operator U2i-i.i that maps this su- 
perposition of two states into a single basis state such 
that 



U2 



1 



(57) 



Since we know that the unknown string z is either in set 
"!v/2-i "-"^ ™ "i^^ wish to measure the outcome of 
the qudit |s)s in the spin basis. We define the projection 
operator [SH 



Ms = \s), {s 
and outcome probability is 
1 



Pr[s] = 



+ 3 



M, 



+ J 



(58) 



(59) 



C. Classical Algorithm 

In this subsection we compare the performance of any 
classical algorithm to the performance of the quantum 
algorithm. 

Claim 3. Any classical deterministic algorithm re- 
quires n(n) oracle queries of the bit positions to solve 
the close Hadamard problem with certainty, even if there 
are no bit errors. A randomized algorithm with bounded 
error probability also requires ^l{n) queries, even if there 
are no bit errors. 

Claim |3] follows from information theoretical consider- 
ations. The goal of the classical strategy is to determine 
which of the A^/2 Hadamard codewords is loaded into the 
oracle. The number of possible solutions is then initially 
n(2"). Whenever a classical strategy performs a query, it 
can eliminate at most half of the remaining possible solu- 
tions, even if there are no errors. To reduce the number of 
possible solutions to a single solution, the classical strat- 
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egy therefore requires at least f2(n) queries"'^. The lower 
bound also holds when the strings loaded into the oracle 
are Hadamard codewords with errors. 

In the next session we discuss how changing the unitary 
operators R and in the algorithm shown in Fig. [1] 
changes the oracle decision problem that can be solved. 



V. ALTERNATIVE ALGORITHM 

The continuously-parameterized, finite-dimensional 
Hilbert space of the spin system inspired an efficient al- 
gorithm for the solution of the close Hadamard problem. 
The group structure of the Hadamard codewords is im- 
plicit in the use of Hadamard operators in the quantum 
algorithm. We now show that this computation model 
can inspire other algorithms. Other unitary operators 
can be employed in the quantum circuit shown in Fig. [T] 
The discrete Fourier transform [sH ] is an obvious alterna- 
tive. We provide a sketch of how the Fourier transform 
changes the group structure of the codewords and point 
to the need for further exploration of problems that could 
benefit from this computational model. 

We replace the operators R and R^ in Fig. [1] with the 
discrete Fourier transform F and F^. The matrix repre- 
sentation of the discrete Fourier transform is expressed 



1 1 

1 iO 



,N-1 



(62) 



V 1 cj(^-i) ••• w(JV-i)(JV-i) J 



where uj = |31j . 

In our analysis of the F-based algorithm, we adopt a 
similar approach to that taken for the H-based algorithm 
and define the 'Fourier codewords' as 



T(^) = log(_i) [V7VF(~) 



(63) 



Similar to the Hadamard codewords, the j*^ Fourier 
codeword is the j"^ row of the matrix T^^^ As an ex- 
ample, we express the N = 8 matrix as 



T(8) ^ 



/ 
1 

Vo -I 





1 

2 

1 

1 

2 



1 

2 

1 

1 

2 





4 

1 

_ 1 

1^ 

2„ 







_ 1 

2 

1 

1 

2 



_ 1 

2 

1 

1 

2 



\ 



(64) 



^ See for example paragraph 6.1 in 
formation theoretic lower bounds. 
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with t]^^ = 01010101. We see that the Fourier code- 
words are not bit strings but rather can be thought of as 
fractional bits. These fractional bits can still be encoded 
into the oracle function U2 given in Eq. ([H]). 

We now define what we term the simple Fourier code- 
word oracle decision problem and show that it can be 
solved in a single query using the modified algorithm. 
Note that we have structured this problem along the 
same lines as the close Hadamard problem with no er- 
rors. 



(JV) 

N/2-1 



and a set of 



Problem 3. Given the string A — 

strings B = Tf \ with k G {Ln \ k ^ N/2 - 1} and a 
string z randomly selected with uniform distribution fi 
such that z C = AU B, the Simple Fourier Code- 
word problem is to determine if z £ A or z E B with the 
fewest oracle queries. 

The action of the algorithm on the input state is ex- 
pressed as 



For z — Tj'^^ and j G 0, 1, . . . , — 1, the output state 
can be shown to be 



(65) 



1*3) 



1 



1 



(66) 



where ^ + j and 
2^2/ 



that 



1 

2 

-s). 



j are modulo s sums in the sense 



We apply 1/2^1 given in Eq. (f57| to the state j^fa) given 
in Eq. (I66p . Measuring the qudit \s)s in the spin basis 
with the measurement operator Ms given in Eq. (j58l) . 
distinguishes whether the encoded string is in set A or set 
B thereby solving the simple Fourier codeword problem 
in a single query. 

The result given by Eq. (j66p is achieved by exploiting 
group properties similar to those expressed in Eqs. (|4T|) 
and for the Hadamard codewords. The columns of 
p(JV) represent the multiplicative cyclic group of order 
N, where the generator is the first non-trivial column 
of F^^). The matrix of codewords T*^^) represents the 
additive cyclic group of order as a result of taking the 
logarithm of F*^^^. 

Each element of the group has the inverse relation 

T,+Tn-,^To, (67) 

and each codeword also obeys the sum relation 



T, + T, 



JV/2-J 



N/2j 



(68) 



where N — j and A^/2 — j are understood to be mod- 
ulo A^ sums. In Fig. [8] we clearly see that, unlike, the 
Hadamard codewords that preserve the superposition of 
two symmetric states, the Fourier codewords preserve the 
superposition of two adjacent states. 
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FIG. 8: For s = 7/2, the Fourier-based algorithm probability 
distributions for the states given by Eq. H66|) for (a) j — 0, 
(b) j — 2 (c) j = 4 and (d) j = 6. The F-based algorithm 
preserves the 'adjacency' of the input superposition, whereas 
the H-based algorithm preserves the 'mirror' symmetry of the 
input superposition as shown in Fig. (6] 



Comparison of the effect of the different operators is 
interesting. The Hadamard codewords preserve symmet- 
ric superpositions, and there are N/2 unique symmetric 
superpositions. The Fourier codewords preserve adjacent 
superpositions, and there are N unique adjacent super- 
positions. 

The exploration of the structure of error cancellations 
along the lines of Eq. ((5(I| using the relationship Eq. ([55)) 
for the Fourier codewords may lead to new problems that 
can be efficiently solved using this model of computation. 
For example, it is natural to apply this error cancellation 
concept to the simple Fourier codeword oracle decision 
problem presented in Problem [3] so that the equivalent of 
Fourier codewords with errors may be included. 



have shown that the optimally-squeezed coherent spin 
state is well approximated by a superposition of two dis- 
crete states. 

Using this idealized state as input into a quantum al- 
gorithm presented in Fig. [1] which is adapted from the 
single-mode continuous variable case, we have discovered 
a new oracle decision problem, which we call the close 
Hadamard problem. This problem is related to the dig- 
ital error correction technique of employing Hadamard 
stings as error-tolerant codewords. 

We have shown that the unrestricted version of this 
problem can be solved in a constant number of queries 
with arbitrarily small error probability, whereas its classi- 
cal counterpart requires fl{n) queries. More interestingly, 
we have also shown the restricted version of the problem 
is solved in a single query with certainty. 

The observed speedup, in the restricted case, results 
from the combination of the group structure of the 
Hadamard codewords, the employment of the Hadamard 
operator H®" as the operator R in Fig. [T] and the use 
of a symmetric superposition of two basis states as the 
algorithm input. The tolerance of errors in this case is 
a direct result of the cancellation of errors that results 
from this combination. 

We have also shown that this algorithm can be gen- 
eralized. We can change the operator R in Fig. [T] to 
the discrete Fourier transform, and loading Fourier code- 
words into the oracle preserves the symmetry of the input 
state superposition in a manner analogous to the close 
Hadamard problem. This pairing between operators and 
the codewords loaded into the oracle offers the promise of 
discovery of new problems that can be efficiently solved 
this way. 

We conclude that the continuously-parameterized rep- 
resentation of quantum dynamical systems having a 
finite-dimensional Hilbert space gives us a model of quan- 
tum computation that inspires efficient solutions of new 
problems. 



VI. CONCLUSIONS 



We have shown that the continuously-parameterized 
coherent spin state representation of a finite-dimensional 
Hilbert space of elementary 1/2 spins gives us a new 
continuous variable model of quantum computation. 
Like continuous variable quantum computation using the 
states of the harmonic oscillator, this spin system is 
amenable to physical preparation with linear rotation 
and quadratic squeezing operators. Unlike the harmonic 
oscillator case, spin squeezing is Heisenberg limited. We 
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